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ABSTRACT 

I investigate the prospects for using the Space Interferometry Mission 
(SIM) to measure the masses of nearby stars from their astrometric 
deflection of more distant sources, as originally suggested by Paczyhski 
and by Miralda-Escude. I derive an analytic expression for the total 
observing time T tot required to measure the masses of a fixed number 
of stars to a given precision. I find that T tot oc r~^ x , where r max is the 
maximum radius to which candidates are searched, or T tot oc ^ in , where 
yU m i n is the minimum proper motion to which candidates are searched. 
I show that Ttot can be reduced by a factor 4 if source availability 
is extended from V s = 17 to V s = 19. Increasing r max and V s and 
decreasing /x min all require a significantly more agressive approach to 
finding candidates. A search for candidates can begin by making use of 
the Luyton proper motion catalog together with the USNO-A2.0 all-sky 
astrometric catalog. However, a thorough search would require the 
all-sky USNO-B proper-motion catalog which is not yet available. The 
follow-up observations necessary to prepare for the mission will become 
more difficult the longer they are delayed because the candidate pairs are 
typically already within 1" and are getting closer. 

Subject headings: astrometry - Galaxy: stellar content - 
gravitational lensing 



- 2 - 



1. Introduction 

Refsdal (1964) pointed out that it should be possible to measure the masses of 
nearby field stars from the astrometric deviation they induce on more distant sources 
as they pass by the latter. To be practical, the two stars must pass within 0(1") of 
each other. Paczyhski (1995,1998) and Miralde-Escude (1996) examined this idea in 
the context of the current rapid improvements in astrometric capability. They made 
rough estimates of the number of mass measurements that could be obtained using 
various ground-based and space-based facilities. Here I re-examine this problem 
specifically guided by the capabilities and requirements of the Space Interferometry 
Mission (SIM). 

The planned SIM launch date is 2005 and the minimum mission lifetime is 5 
years. In order to carry out mass measurements, two steps must be completed prior 
to launch. First, one must identifify candidate pairs of stars from a proper- motion 
catalog: a nearby "lens" star must be found that is likely to pass sufficiently close 
to a more distant "source" star to cause a large deflection of light and so permit 
a precise measurement of this deflection. Second, given the quality of the catalogs 
that will be available in the near future, it will generally not be possible to predict 
which of the candidates will be the best to make precise mass measurements 
with a modest amount of observing time. Rather, it will be necessary to perform 
follow-up observations of these candidates prior to the event in order to determine 
the impact parameter (angular separation (3 at the point of closest approach). 
Typically, the candidates are already closer than 1", often much closer. Moreover, in 
many cases, one star is substantially brighter than the other. Hence, the follow-up 
observations will usually require adaptive optics or the Hubble Space Telescope. 
These requirements will grow more severe as time passes. In brief, preparation for 
mass measurements using SIM should proceed without delay. 

Because SIM observing time comes at a high premium, my approach is to 
rank candidates by the amount of observing time that is required to make a mass 
measurement of fixed precision. I then use this framework to characterize and 
evaluate various selection strategies. 

The probability that it is possible to measure the mass of a given foreground 
star grows monotonically with its proper motion and is linear in the proper motion 
in most cases. Hence, a survey based on an ideal star catalog (not affected by 
magnitude limits or crowding) would investigate foreground stars down to some 
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minimum proper motion /i m i n - On the other hand, for stars of sufficiently low 
luminosity, the magnitude limits of the underlying catalog will impose an effective 
distance limit, r max . Thus, it is important to consider both forms of selection. In 
practice, the actual selection process may also be affected by crowding, but I will 
not consider crowding explicitly in this paper. Rather, one may think of crowding as 
imposing an indirect constraint on yU min or r max . 

I derive simple expressions for the total observing time T to t needed to make 
M mass measurements. For fixed J\f, I show that T tot oc r~^ x for distance limited 
surveys, and T tot oc /i^ nin for proper-motion limited surveys. Hence, minimization 
of the observing time requires pushing r max out as far as possible or pushing /i min 
as low as possible. The sample of candidates will then have on average smaller 
proper motions meaning that they are even closer on the sky today, thus making 
follow-up observations even more difficult. In addition, I show that by extending 
the available sources from V s = 17 to V s = 19, one can decrease T tot by a factor of 
4 despite the lower flux from these fainter sources. However, to determine which 
of these fainter sources are really usable requires a much more precise estimate of 
their expected impact parameters, that is, even more precise measurements of their 
current positions despite the larger disparity in the source/lens flux ratio. Obviously 
these measurements will also become more difficult with time. 

To carry out a search for candidates it would be best to begin with an all-sky 
proper motion catalog. Such a catalog is currently being prepared by the US Naval 
Observatory (D. Monet 1999, private communication) but has not yet been released. 
In the meantime, one can make a good start using the Luyton (1979) proper motion 
catalog in combination with the USNO-A2.0 astrometric catalog (Monet 1998). I 
briefly describe how to carry out such a search. 

2. Required Observation Time For an Individual Lens 

Consider a nearby star ("the lens") of mass M and at distance r that 
passes within an angle (3 of a more distant star ("the source") at r s . The 
source light will then be deflected by an angle a = 4GM/(/3rc 2 ) at the 
point of closest approach. Consequently, the source will appear displaced by 
a = a(l — r/r s ) = 4GM7r re i/(AU ( 5c 2 ) relative to the position expected in the absence 
of lensing. One can therefore determine the mass of the lens by measuring this 
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displacement, provided that (3 and the relative parallax 7r rc i are known. Note that 

M f Br \- 1 
a = 80uas— — -— . 1 
P M Q V100AU/ v ; 

Assuming photon-limited astrometry measurements, the total amount of 
observing time r required to achieve a fixed fractional error in the mass measurement 
then depends on three factors. First, the measurement is easier the bigger a: the 
factional error for fixed observing time falls as oT x and so for fixed fractional error, 
r oc oT 2 . Second, the measurement is easier the brighter the source magnitude V s , 
t oc 10 0AVs . Third, the time required depends on the geometry of the encounter. 
The geometry can be described in terms of the angular coordinates A) of the 
source-lens separation vector at the time (t = 0) of the midpoint of the mission and 
the angular displacement fiAt of the lens relative to the source during the course 
of the mission. Here (3 is the source-lens separation at the time of closest approach 
(t — to), fj, is the relative source-lens proper motion, At is the duration of the 
mission, and A = —//to- I therefore write, 

where 7 is a function to be described below, and where «o and To are convenient 
normalization factors. For definiteness, I will take the required mass precission to 
be <7m/M = 1% and will arbitrarily adopt «o — 100 //as. I will normalize 7 so that 
it is effectively the number of equal-duration measurements that must be made. I 
characterize SIM astrometry as requiring 1 minute to achieve 40 //as precision in 1 
dimension at V s — 17. Then 

To estimate 7, I consider sets of observations over the angular interval [A_, A+] 
where A± = A ± //At/2, and solve simulteously for six source parameters: the 
two-dimensional angular position at the midpoint, the two-dimensional proper 
motion, the parallax, and a. Even though very little information can be obtained 
about a from astrometry measurements parallel to its direction of motion, I include 
such measurements in order to be sensitive to other kinds of apparent source 
acceleration (e.g. gravitational). Without such a check, the mass measurement could 
not be considered reliable. I then optimize these observations for the measurement 
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of a. Generally, the optimum configuration has roughly equal total exposure times 
at the point of closest approach and near the beginning and end of the experiment, 
and has no observations at other times. I take At to be 5 years, but the results 
would be the same for any value provided that At 3> 1 yr. I find that 7 achieves a 
minimum, 7 ~ 10 when A_ < —2/3 and A + > +2(3. For example, 7(0,2;) = 10 for 
x > 4. Some other indicative values are 7(0, 2) = 21, 7(0, 1) = 99, 7(0.75, 2.5) = 24, 
7(0.25, 1.5) = 39. Thus, if the source does not move by a distance at least equal to 
the impact parameter on either side of the lens during the course of the observations, 
7 becomes very high. In my analysis below, I will incorporate the exact values of 7 
for each configuration. But qualitatively one can think of 7 as being 

where is a step function, and 7* = 10. 



3. Observing-time Distribution 

From the previous section, a star with M = M & , r = 100 pc, (3 = 1", and 
a minimal 7 would require about 420 hours of observation time for a 1% mass 
measurement. Hence, it is prudent to consider how one might find pairs of stars 
with the most favorable characteristics. I begin by writing down the observing-time 
distribution for an arbitrarily selected sample but with lenses of fixed mass M, 

^ = J d 3 rd 3 r s dV s dv ± n(r)n s (V s ,r s )f(r,v ± )S(r,r s ,v ± ,V s ,...)5p 

(5) 

Here n(r) is the number density of lenses as a function of their position, n s (V s ,r s ) 
is the number density of sources as a function of their magnitude and position, 
v± = rfi is the transverse speed of the lens relative to the observer-source line of 
sight, /(r, v±) is the transverse speed distribution as a function of position, S is 
the selection function (with possibly many variables in addition to those explicitly 
shown), b = rf3, £ = rX, and 6 is a Dirac delta function. Equation (0) cannot be 
simplified without additional assumptions. As I introduce these assumptions, I will 
briefly outline their impact. Some of the simplifications will then be discussed in 
greater detail below. 

I first assume that r<r s . This is an excellent approximation for disk lenses 
although it is not as good for halo lenses. It has two simplifying effects: a — > a, 
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dv ± dV s J df3d\drr 2 nQ(r max -r)5(T - r) J dQ, (j)(V s , tyf(v ± , r). (7) 



so r = r(Vs, M, 6, £, v±At), and the 3-space density of sources n s can be replaced 
with the projected surface density <p(V s ,Q), where Q is position on the sky. (To be 
more precise <ft is the density of sources in the neighborhood of the lens postion Q.) 
Second, I assume that the product of the selection function and number density can 
be written, 

n(r)S(r, r s , v±, V s , ...) -> n6(r max - r), (6) 

where n is now assumed to be uniform, and r max is a maximum search radius. In 
fact, this is an oversimplification. A major focus of the present study is to determine 
what effect the selection function has on the observing-time distribution. The best 
way to do this is to begin with this simplified picture. With these assumptions, 
equation (|5|) can be written, 

dM 
~dT " 

In this form, the integration still cannot be factored because of the correlation 
between the speed distribution of the lenses and the distribution of sources. I 
therefore assume that f(v±,r) — > f(v±), i.e., that the speed distribution does not 
depend on position. This is actually a very minor assumption, provided that the 
speed distribution is taken to be the average over the Galactic plane where the 
majority of the source stars are. Two of the integrals can then be evaluated directly, 
and equation (0) becomes, 

dU 
If 

where </>(V^ s ) = / dQ(p(V s , Q) is the luminosity function integrated over the entire sky. 
Equation (|§p already contains an important result: the number of lenses available for 
measurement at fixed observing time is directly proportional to r max , the physical 
depth to which they are searched. 

To further evaluate the integral, first define, 

G(i;x) = JdyStf-ifax)]. (9) 

Then the integral can be written 
dj\f 



nr max / dbd£dV s dv ± f{v ± )<f>{V s )8[T-T{V s ,M,b,e,v ± At)}, 



^r max J dVdV s (j>{V s )8[T - t(V s , T)]H(T), (10) 
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3.1. Analytic Estimate 

I will evaluate equation (|It]) explicitly in § 3.2 below. However, it is also 
instructive to make an analytic estimate of this equation with the help of a few 
approximations. First, I assume that all the sources have the same magnitude, 
4>(V S ) = N8(V S — 17), where N is the total number of source stars. Hence, 



dT T \T 0/ 
Second, I use the approximation (§) to estimate G, 

G( T ,x) = (x-2)e(x-2)8(j*-j), (7* = 10). (13) 

Third, I take f(v±) = 5(v± — t>*), where v* is a typical transverse speed for the lens 
population. Then 



H(T) = T-^bo^ 



, T \ 1 /2 1 2GM 

v*At - Ab 

' 7: 



&o = r , (14) 



where I have suppressed the function that limits the range of validity to 
T < •y*(aoC 2 v*At/8GM) 2 . Combining equations (|T2| ) and (p^|), I obtain, 



cW 2GpNr^v*At 



«oc 2 v*At\ 2 



(15) 



dT ( 7 ,TT ) 1 /2 aoC 2 ' 
where p = nM. The limiting condition in equation (|T5| ) comes from assuming that 



the first term in brackets in equation (14) is much greater than the second. Equation 



(|i~5|) tells us that the observing-time distribution depends on the type of lens only 
through its mass density p, its typical velocity t>*, and the cutoff which scales as 
{vjMf. 

A sensible observing strategy will naturally focus on the lenses that require the 
least observing time. I therefore consider a program that measures the masses of all 
lenses requiring observing times less than some maximum, T max . The total observing 
time T tot can then be expressed as a function of the total number of stars observed, 
A/", and of the other parameters: 

T M =[ d TT§, M = l ,T§, (16) 



1 . rZ (AGpNr m!at y*At 



-2 



3 V a c 2 1 
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Equation ([17]) is one of the major results of this paper. It states that the total 
observing time required to measure the masses of a fixed number of lenses scales 
inversely as the square of the search radius r max of the sample. Given the premium 
on SIM time, this result implies that the search should be pushed to as large a radius 
as possible. I discuss the prospects for doing this in § 4. 

The total time can be written out explicitly 

where I have assumed a mission lifetime of At = 5 yrs and normalized the transverse 
speed to a typical disk value v* = SSkms" 1 and the density to approximately 1/3 of 
the local stellar disk density (Gould, Bahcall, & Flynn 1997). That is I consider that 
one is interested in one (or perhaps several) subsets of the whole disk population. I 
have also assumed a total of N = 10 8 stars at V s = 17 over the whole sky (Mihalas 
& Binney 1981). This estimate incoporates a maximum observing time per object, 

3 

-^rnax ^y--^~tot > (1^) 

which must be well under the cutoff in equation flT5| ) given by 

/ v*At \ 2 ( M \- 2 , , 

^ = 13 hours (— ) (— ) . (20) 

If T cut ^ T max , then the scaling relation ( |I~7D [M oc T^) is no longer satisfied. See 
Figure [l|, below. The cutoff is satisfied for low mass disk stars (assuming only 5 mass 
measurements are desired) but becomes more difficult for higher masses. 

Another important feature of equation (|l~7|) is that T tot oc 10 0A( - Va ~ 17 ^T / N 2 . 
Thus, if we compare V s = 17 and V s = 18, the latter are 2.5 times fainter and so 
require 2.5 times greater 10°' 4 ^ s_17 ^To, the observing time for a single astrometric 
measurement of precision oq. On the other hand, there are approximately 1.9 times 
as many stars (Milhalas & Binney 1981) and so T tot is actually smaller by a factor 
~ 0.7. It should be noted, however, that the shorter observing time comes about 
because the impact parameter b is typically 1.9 times smaller. In § 4, I will discuss 
the prospects for recognizing when such close encounters will occur. 
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3.2. 



Numerical Estimates 



To test the estimates derived in § 3.1, I continue to approximate as a 
S function, but otherwise carry out the full integration indicated by equations 
(|10|) and flTTl). I take the velocity distribution to be a two-dimensional Gaussian 
with (one-dimensional) dispersion typical of foreground objects in the Galactic 
plane: a 2 = crfj/A + cry/A + a 2 v /2, where a v = 34kms _1 , ay = 28kms _1 , 
and aw = 20kms _1 . This yields a = 26kms~ 1 for which the mean speed is 
f* = (7r/2) 1 / 2 cr = 33kms _1 (as used in § 3.1). Figure [I] shows the results for 
M = M (bold curve) and M = 0.1 M Q (solid curve). The agreement with the 
analytic prediction from equation (0) (dashed line) is excellent. Equations ( |T9"D and 
(|2g) predict that the cutoff should be at TV ~ 1.5(M/M ) _1 , or at TV ~ 1.5 and 
TV ~ 15 for the two cases shown. In fact the actual values are about 2.5 times higher. 
Most of this difference (a factor of 2) is due to the fact that the velocity distribution 
is not a 5 function, and the higher-speed stars are more likely to be candidates and 
are less affected by the threshold. 

Figure |2| shows the same quantities for six different magnitude bins of the 
luminosity function (Mihalas & Binney 1981), and M = M & . The V s = 17 curve 
(same as in Fig. [l|) is shown as bold dashed line, and the others V s = 14, 15, 16, 18, 19 
are shown as solid lines. The curves can be separately identified by noting that the 
cutoff increases with magnitude. The upper bold line shows the result of combining 
all of these while the lower bold line shows the result of combining the four bins 
with V s < 17 together. Note that each of three bins V s = 17, 18, 19 contribute about 
equally (for T tot ^ 100 hours). This is because the longer integration times required 
for the fainter sources are compensated by the fact that they are more numerous 
and hence closer on average to the lenses. 



As I discussed in the introduction, in some regimes the selection function is best 
described as a cut on distance and in other it is best described as a cut on proper 
motion. So far, I have focused on selection by distance. See equation @. Had I 
instead selected on proper motion, 



3.3. Proper-motion selection 




(21) 



-10- 



then equations (|H]) and (|TTD would be replaced by 
dj\f 



n J drdTdV s <P(V s )5[T - r(V s , T)]H{T; r// min ), (22) 



dT 
where 

H(T;u) = £dv ± f(v ± ) J djdbbG^r,^-)*^ - (^)"V 

Carrying through the derivation, one obtains the analog of equation (|1 

l^ r3 MGpN(vl)A t V \ 
3 V /iminaoc 2 



(23) 



T tot = ^(^™^) 7 . Toi (24 ) 



where (f^) is the mean square transverse speed. That is, equations ( |17|) and (|4]) are 
identical except rmaxW* — > (vj_) / fi m i n . For a Gaussian, (t> 2 ) = 2a 2 , so this relation 
can be written 

1 r 90 pc — ^ _ 1 Y\ (25) 



7r /i min 33kms _1 VlOOmasyr 

I find that with this substitution, the curves produced by equations (|22"1) and (|23| ) 



are almost identical to those produced by equations (|T0|) and (|Tl| ) except that the 
cutoffs are increased by a factor 1.2. Thus, proper-motion selection and distance 
selection produce essentially the same results, provided they are converted using 
equation (|25|). 



Identification of Pairs 



From equation ([18]), the total observing time required to measure the mass of 
a fixed number of stars declines as r~ 2 x . From Figure || one sees that including the 
magnitude bins V s = 18, 19 is roughly equivalent to increasing the total number of 
sources iV by a factor 2, and from equation (|T8|), the observing time is reduced as 
iV~ 2 ~ 0.25. This estimate is confirmed by the offset between the two bold curves in 
Figure 0. 

Hence, if it were possible to push out to fainter sources and more distant lenses, 
it would certainly be profitable to do so. I therefore now investigate the constraints 
governing the identification of lens-source pairs. The basic problem is that if these 
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pairs are to be close enough for astrometric microlensing in, say, 2008, then they are 
already very close. That is, their separation A# now is 



Thus, it would be difficult to conduct a large-scale survey for such pairs. 

Fortunately, the US Naval Observatory is soon planning to release an all-sky 
proper motion survey, USNO-B (D. Monet 1999, private communication). Even so, 
the identification of pairs is not trivial, and becomes more difficult both for fainter 
sources and larger lens distances (or lower proper motions). 



I begin by considering the problem of candidate pair identification when 
the proper- motion survey data conform to "typical" specifications. In fact, the 
underlying data sets are heterogeneous, with substantially longer baselines in the 
north than the south. For definiteness, I will consider proper motions derived 
from 2 epochs, one in 1955 and the other in 1990. This baseline is appropriate 
for the northern hemisphere. The anticipated proper-motion error is 4masyr~ 1 , 
corresponding to 100 mas errors in each position measurement. This implies an error 
of about 120 mas in the predicted positions of the source and the lens in 2008, or 170 
mas error in their relative position. (Generally, only the error in one direction - that 
of the impact parameter - comes into play.) Is this good enough? Let us suppose 
that all pairs requiring T < T max = 200 hours are to be observed. From equation (|I9"D 
this corresponds to J\f ~ 6. For 7 = 10, equation (0) implies a < 115 //as lo - 2 ^ -17 ), 
and so from equation (TO) 



Hence, for M ~ M Q , /3 max is greater than 170 mas even for r = 200 pc or V s = 18. 
Thus while it would still be necessary to do additional astrometry in order to prepare 
for the observations, relatively few candidates would be rejected by this astrometry. 
By contrast, for M = 0.2 M , r = 200 pc, and V s = 18, /3 max ~ 30 mas. In this case, 
it would be necessary to examine about 6 candidates drawn from the proper-motion 
survey to find one suitable for a mass measurement. 




(26) 



4.1. Unblemished Survey Data 




(27) 
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In the southern hemisphere, the baselines are shorter, so the proper-motion 
errors are about twice as big. Hence, about twice as many candidates need to be 
examined. 

4.2. Compromised second epoch 

If the time of closest approach is 2008, then at the time of the second epoch 
of the proper motion survey, say 1990, a typical foreground star will be separated 
from the line of sight to the source by of order 100 AU. This corresponds to 1" at 
r = 100 pc. As discussed in § 3.1 in the analysis of Figure ^ lens candidates tend to 
be moving faster than the population as a whole, so the actual typical separation 
will be somewhat larger. Nevertheless, these values are close to the resolution limit 
of the surveys and become even less favorable at greater distances. In addition, 
bright lens stars will entirely blot out a substantial region around them on the survey 
plates, preventing the detection of candidate source stars at all. For example, I find 
that on the Palomer Observatory Sky Survey (POSS), V = 8 stars (the approximate 
completeness limit of the Hipparcos catalog) tend to blot out a region with a 10" 
radius. 

However, even the complete loss of the second-epoch positions of candidate 
sources is not crippling. The proper motion of the bright lens candidate can still be 
measured, and its position in 2008 predicted. This region can then be examined on 
the first epoch plates for potential candidates (assuming that the lens is not bright 
enough to have blotted out this region even at this earlier epoch). Of course, in 
the intervening ~ 50 years, these candidates will have moved, but probably not by 
very much. For example, at high latitudes (|6| ^ 20°), disk sources typically lie at 3 
disk scale heights or about | csc6| kpc. Hence, in 50 years, they will typically move 
300 1 sin &| mas which even at b = 90° is not much more than the error in the expected 
position (see § 4.1). Closer to the plane, the motion will be even less. Thus, without 
a second epoch, more candidates will have to be examined at high latitudes (but 
these contain a minority of the candidates anyway) and there will be hardly any 
effect at low latitudes. 



-13- 



4.3. What can be done now? 

The USNO-B catalog has not yet been released. However, it is still possible to 
begin the search for candidates using the Luyton (1979) proper-motion catalog in 
combination with the USNO-A2.0 all sky astrometric catalog (Monet 1998). The 
conditions of such a search are fairly well described by § 4.2. 

Based on a comparison to Hipparcos data, I.N. Reid (1999 private 
communication) estimates that the Luyton (1979) proper motions are typically 
accurate to about 10 mas yr _1 , and that the positions are accurate to a few 
arcseconds. By identifying the Luyton (1979) stars in the USNO-A2.0 survey, 
one could fix the ~ 1955 positions to ~ 100 mas, and hence the 2008 positions to 
~ 500 mas. One could then search the USNO-A2.0 catalog for background stars 
whose ~ 1955 positions lay along the 2008 path of the Luyton star. As discussed in 
§ 4.2, these stars could be expected to move about 300 1 sin &| mas which is generally 
small compared to the uncertainty in the position of the foreground star. Hence, 
2008 source-lens separations could be predicted to ~ 0."5. For pairs that were 
sufficiently close, the separation could be measured on the POSS II plates to refine 
the prediction. Follow-up observation could then be made of those pairs surviving 
this test. 

5. Stellar Halo Lenses 

Halo stars are about 500 times less common than disk stars (Gould, Flynn, & 
Bahcall 1998), i.e., p ~ 6 x lO _5 M pc -3 , and they are typically moving about 5 
times faster. Let us suppose that they could be spotted to r max = 1 kpc (see below), 
and let us take N = A x 10 8 in accord with the discussion of Figure 0. Then, from 
equation (0), it would be possible to obtain the masses of 5 halo stars with about 
150 hours of observation. 

At r max = 1 kpc, it is still appropriate to approximate the density of the stellar 
halo as uniform. However, it is no longer appropriate to treat the sources as being 
infinitely far away. As mentioned in § 4.2, at b = 90°, typical disk sources are at 1 
kpc. However, the disk sources are farther away at lower latitudes where there are 
the greatest fraction of sources in any case. Hence, given the level of approximation 
of the present study, I will ignore this modest correction. 
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Of course, the radius within 1 kpc contains of order 300 times more stars than 
the radius within 100 pc, so identifying a relatively complete sample of halo stars 
seems like a formidable task at first sight. In fact, for stars of fixed color (and 
approximated as black bodies - or at least as deviating from black bodies in similar 
ways), we have 

t croBB = k {B - R)o 10- - 2Ro t i- 1 , (28) 

where km-R) is a constant that depends on the dereddened color, Rq is the 
dereddened magnitude, and t cross is the time it takes the star to cross its own radius. 
For halo stars, t cross ~ 10 3 s, which is significantly different from the value for other 
common classes, 2 x 10 2 s for disk white dwarfs, 5 x 10 3 s for thick-disk stars, and 
2 x 10 4 s for main sequence stars. Thus it should not be difficult to find halo star 
candidates in a proper-motion catalog with colors. The sample will be somewhat 
contaminated with fast-moving thick-disk stars, but these are also of considerable 
interest because of their low metallicity. 

A more significant problem is that if the survey is limited to V ~ 20, then at 1 
kpc, the bottom of the luminosity function My > 10 is not detectable. These fainter 
stars contain about half the spheroid mass, implying that the above estimate of the 
observation time required to measure 5 masses should be multiplied by 4 to about 
600 hours. Of course, if one were willing to settle for 3% measurements in place of 
1%, this estimate would come down by an order of magnitude. 



(29) 
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Fig. 1. — Total number of mass measurements M as a function of total required 
observing time T tot . The solid curve is for M = 0.1 M & and the bold curve is for 
M = M Q . The dashed line is the analytic approximation (but without the cutoff) 
given by eq. (jlB"!). 



Fig. 2. — Total number of mass measurements M as a function of total required 
observing time T tot , for mass M = M & . The V s = 17 curve (same as in Fig. p]) 
is shown as a bold dashed line, and the others V s = 14, 15, 16, 18, 19 are shown as 
solid lines. The curves can be separately identified by noting that the cutoff increases 
with magnitude. The upper bold line shows the result of combining all of these bins 
(14 < V s < 19) while the lower bold line shows the result of combining the four 
brightest bins (14 < V s < 17). 



